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We study thermal radiation from nonlinear (x^) photonic cavities coupled to external 
channels and subject to incident monochromatic light. Our work extends related work on 
nonlinear mechanical oscillators [Phys. Rev. Lett. 97, 110602 (2006)] to the problem of 
thermal radiation, demonstrating that bistability can enhance thermal radiation by orders 
of magnitude and result in strong lineshape alternations, including “super-narrow spectral 
peaks” ocurring at the onset of kinetic phase transitions. We show that when the cavi¬ 
ties are designed so as to have perfect linear absorptivity (rate matching), such thermally 
activated transitions can be exploited to dramatically tune the output power and radiative 
properties of the cavity, leading to a kind of Kerr-mediated thermo-optic effect. Finally, 
we demonstrate that in certain parameter regimes, the output radiation exhibits Stokes and 
anti-Stokes side peaks whose relative magnitudes can be altered by tuning the internal tem¬ 
perature of the cavity relative to its surroundings, a consequence of strong correlations and 
interference between the emitted and reflected radiation. 
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Driven nonlinear oscillators, including optical,- optomechanical,-- and MEMS^ resonators, 
have been studied for decades and exploited for many applications, from mass detection^ to sens¬ 
ing^ and tunable filtering. 7 ' When driven to a non-equilibrium state, these systems can exhibit a 
wide range of unusual thermal phenomena,- leading for instance to cooling and amplification of 
thermal fluctuations in optomechanical systems,^ generation of squeezed states of light in Kerr 
media,- and stochastic resonances.— Previous studies of Duffing oscillators have also identified 
novel effects arising from the nonlinear interaction of coherent pumps with thermal noise,— 
leading to phase transitions and lineshape alterations that were recently observed in a handful of 
systems, e.g. mechanical oscillator s 1 - ^— and Josephson junctions.-— 

In this letter, we study thermal radiation effects in optically driven x' 31 photonic cavities cou¬ 
pled to external channels. We demonstrate that in certain parameter regimes, bistabilityi2dJ in 
photonic resonators leads to thermally activated transitions that amplify thermal fluctuations by or¬ 
ders of magnitude and cause dramatic changes in the cavity spectrum, analogous to noise-induced 
switching in mechanical oscillators.— We find that when the photonic cavity is critically coupled 
to the radiation channel (enforced by designing the cavity to have equal dissipation and radiation 
rates),— the coherent part of the output power varies dramatically with temperature, leading to a 
kind of Kerr-mediated thermo-optic effect. A simple perturbative analysis also shows that outside 
of the bistability region, the interaction of the coherent drive with thermal noise leads to amplified, 
Raman-type Stokes and anti-Stokes side peaks in the radiation spectrum, the relative amplitudes 
of which depend on a sensitive interference between the externally incident and reflected thermal 
radiation. Related phenomena have been long studied^ and more recently observed^^ in the 
context of driven nonlinear mechanical oscillators as well as resonators based on rf-driven Joseph¬ 
son junctions-^, and microscopic theories have also been used to describe related optical effects 
in the quantum regime.— Our work is an extension of these studies to the particular problem of 
thermal radiation from photonic resonators. As we show below, additional considerations arising 
in the case of radiation from cavities but absent in mechanical oscillators or bulk media, such as 
the strong coupling of the cavity to an external channel, dramatically impact the outgoing radia¬ 
tion. The ability to tune the radiation properties of resonators via temperature and optical signals 
offers potentially new avenues for applications in nano-scale heat regulation,— detection,— recti¬ 
fication,photovoltaics,— or incoherent sources.— We propose a practical, photonic structure 
where these effects can arise near room temperature and at mW powers. 

The system under consideration belongs to the class of nonlinear photonic resonators depicted 


2 




FIG. 1. Schematic of a wavelength-scale silicon ring resonator of radius R = 4.4//m, height h = 220nm, 
and width w = 350nm, coupled to a silicon waveguide (channel), both on a silica substrate. Also shown is 
the E y mode profile of a resonance designed to have azithmutal number m = 25, wavelength A = 1.5pm, 
radiative lifetimes > 10 6 , and relatively large nonlinear coupling coefficient a = 0.032 X ( 3 W(8e 0 A 3 ). 
The loss 7 d and waveguide-coupling 7 e rates are much larger than the corresponding radiation rate. 


in Fig. |H involving a cavity coupled to an external channel (e.g. a waveguide). The description 
of thermal radiation in this system can be carried out via the coupled-mode theory framework—, 
which we recently employed to study thermal radiation in a related system^ but now extend to 
consider the addition of a coherent pump. The equations describing the cavity mode a are given 
by:— 


^ = [i(w 0 - a|a| 2 ) - 

s- = -s+ + 


7]a + y/^Td&i + v / 2% s -h 


( 1 ) 

( 2 ) 


where |a | 2 is the energy of the cavity mode and |s ±| 2 are the incident (+) and outgoing (—) power 
from and into the external channel, respectively. The latter arises due to dissipative noise inside 
the cavity £ d as well as externally incident light s + consisting of both thermal radiation £ e and 
a monochromatic pump s p exp(icu p t). The dynamics of the cavity field are described by its res¬ 
onance frequency uq and decay rate 7 = 7 e + Jd, which includes linear absorption 7 d as well 
as decay into the external channel 7 e . The real and imaginary parts of the nonlinear coefficient 
a = |cc 0 / e 0 x (3) \£\ 4 /(f s\E\ 2 ) 2 depend on a complicated overlap integral of the linear cavity 
fields,— and lead to self-phase modulation (SPM) and two-photon absorption (TPA), respectively. 
We mainly focus on the effects of SPM (real a > 0) since we find that TPA leads to thermal broad¬ 


ening of the kind explored in Ref. 28J. Both internal and external thermal sources are represented 


by stochastic, delta-correlated white-noise sources £ e and £ d satisfying (assuming 7 u; 0 ), 




(3) 
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where (...) denotes a thermodynamic or ensemble average, and ©(a;, T) = Hcu/ _ i) i s 

the mean energy of a Planck oscillator— at local temperature T; the temperatures of the internal 
and external baths are denoted as T d and T e , respectively. Above, we assumed fn/k n -C T e , T d al¬ 
lowing us to ignore the frequency-dispersion and temporal correlations (colored noise) associated 
with 0. (Note that in the limit hwo/k^T — > 0 one obtains the classical result 0 —>• k B T.) 

Thermal amplification and power tunability .— We show that bistability can amplify thermal 
fluctuations and lead to enhanced, temperature-tunable emission from the cavity. We begin by 
reviewing a number of key features of the system in the absence of thermal noise, whose contri¬ 
butions are considered perturbatively due to the generally weak nature of thermal noise, |s p | 2 A> 
pk]]T. The steady-state cavity field ao due to the pump is given by the well-known cubic equa¬ 
tion 


A 


a|ao| 

7 


+ 1 


a|ao| 

7 


= 2C, 


(4) 


where ( = a|s p | 2 7 e / 7 3 is the effective nonlinear coupling associated with the pump and A = 
is the dimensionless detuning. Equation © describes a number of extensively studied non¬ 
linear effects,—^- including bistability arising in the regime A < —\/3 and (A 1 ) < ( < (A 2 ), as 
illustrated by the hysterisis plot on the inset of Fig. ©a) which shows the dimensionless cavity 
energy x = a|a | 2 /7 as a function of (. 

An effect that seems little explored but that plays an important role on the thermal properties 
of this system is perfect absorption, which occurs when a photonic cavity is driven on resonance 
and its dissipation and radiation rates are equal, also known as rate matching.— In the presence of 
nonlinearities, the cavity frequency and hence the absorbed power depend on For instance, in 
the non-bistable regime, the output power varies slowly with as illustrated by the green curve in 
Fig. ©a) for A = —1, increasing and then decreasing as ( —>■ | A/21, at which point the cavity and 
pump frequencies are in resonance, i.e. a|a 0 | 2 /7 = —A. Bistability can lead to a more pronounced 
dependence on (: the two stable steady states experience different frequency shifts and hence loss 
rates, and ultimately which state is excited in the steady state depends on the specific initial (or 
excitation) conditions.—Figure ©a) shows the steady-state output power |s_| 2 as ( is adiabatically 
increased (solid lines) from zero and above the critical point <AA, for multiple A. The dashed blue 
line shows the power as ( is adiabatically decreased below ( (2> for the particular case A = —2.5, 
demonstrating that only the upper branch experiences perfect absorption, occuring at = A/2 
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FIG. 2. (a) Output power |s_| 2 normalized by the input power |s+| 2 of the pumped system described in 
Fig. SI in the absence of thermal noise and as a function of £, for different values of detuning A = 

The top left inset shows a hysterisis plot of the energy a|ao| 2 /7 as a function of £, the solutions of ©, for 
the particular choice of A = —2.5 while the bottom inset shows the corresponding potential energy U as a 
function of the cavity energy for two different £ = | A/2| and Q 2 '. (b) The same normalized output power 
|s_| 2 /|s + | 2 as a function of temperature (tQ 2 , where Qr = aO(ujQ. Tjy'e/y 2 and Q = cjo / t , for different 
values of A and £ < Q i2) . The insets illustrate the change in the coherent (left) and thermal radiation (right) 
spectra. Both internal and external baths have equal temperatures Xd = T e = T. 


and marked by the white circle. The corresponding change in the output power as the system 
transitions from the lower (A) to the higher (B) energy state at (/T i s given approximately by: 

I 12 A 7e-7d \ (A + xi) 2 — (A + X 2) 2 

1+1 V 7e + 7d/ [1 + (A + 2 t) 2 ][1 + (A + x 2 ) 2 ] K) 

where x\ = — |(2A + \/A 2 — 3) and X2 = — 2(A + xi) are the cavity energies associated with the 

lower and higher energy state, respectively. Given ©, one can show that the difference in output 
power is largest under the rate matching condition 7 e = 7 d and at A ^ —7/3, decreasing with 
smaller or larger detuning. 

The presence of noise complicates this picture due to finite-temperature fluctuations which 
cause the system to undergo transitions between the two states, where the rates of forward/backward 
transitions are a complicated function of the potential energy U and temperature of the systcm.-AA 
(For convenience and without loss of generality, we take both thermal baths to have the same tem¬ 
perature T.) In particular, U is obtained by integrating the steady-state equation associated with 
the cavity energy ^ = [(A + x) 2 + l]x — 2£ = 0 with respect to x = a|a| 2 / 7 - Examples of 
U are shown on the lower inset of Fig. [2 a) for two values of £. Thermally activated hopping 
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leads to significant enhancement of amplitude fluctuations, which manifest as large changes in the 
radiation spectrum of the cavity. This is illustrated by the top inset of Fig. Ob), which shows the 
thermal spectrum of the power |5s_(cu)| 2 for the particular choice of A = —2.5 and ( = 1.3 and 
for multiple values of CtQ 2 , where for convenience (below) we have introduced the dimensionless 
effective thermal coupling ( T = a0(w o , T) 7 e / 7 2 and cavity-lifetime 0 = cu 0 /7-— Such enhance¬ 
ments were predicted to occur and recently observed in nonlinear mechanical oscillators,— 14 
where the authors showed that at special the system undergoes a so-called kinetic phase tran¬ 
sition associated with equal rates of forward/backward hopping and exhibits a “supemarrow” and 
highly amplified spectral peak. Interestingly, we find that in the case of optical resonators, thermal 
amplification can be accompanied by a significant decrease in the coherent output power despite 
the fact that |s p | 2 xk]>T, a consequence of perfect absorption. In particular, operating under 
rate matching and near Q i2} allows for temperature to initiate transitions (A) = (B), leading to 
significant changes in |s_| 2 with respect to T. Essentially, as ( —> </ 2 / the potential barrier 
separating the lower x\ from the higher x 2 energy states begins to dissapear, resulting in increased 
rate of forward transitions and hence larger absorption. 

These features are illustrated in Fig. |2jb) which shows the total output power as a function 
of CtQ 2 for different combinations of (' and A. (We found numerically that for a given ( and A, 
changing either or Q 2 while leaving CtQ 2 unchanged leaves |s_| 2 /|s + | 2 unaltered.) Noticeably, 
while the change in the output power is gradual in the non-bistable regime (A > —a/3), there is a 
significantly stronger dependence in the bistable regime—the slope becomes increasingly sharper 
as ( —> </ 2 ) and ( T —> 0 since it becomes increasingly easier for lower T fluctuations to induce 
hopping unto the higher-energy state. At sufficiently large (t, |s -| 2 is found to increase with 
increasing (t as the cavity field no longer probes the hysterisis regime. While the maximum 
change in |s_| 2 can be estimated from the steady-state analysis in the absence of noise (with the 
largest change occuring for A ~ —7/3), its dependence on ( T is a complicated function of (' and 
A. For instance, for A = —2.5 and ( = 1.4 [blue line in Fig. Ob)], the sharp decrease in output 
power occurs at ( T Q 2 ~ 10 and yields a slope ]~p- ~ 0.05. 

Side peaks. — We now show that the radiation spectrum also exhibits other interesting features, 
including the emergence of Raman-type Stokes and anti-Stokes side peaks previously observed in 
driven mechanical oscillators Interestingly, we find that in our photonic resonator, the pres¬ 
ence of the external channel dramatically alters the relative amplitudes of the side peaks, e.g. 
leading to a symmetric spectrum when the two baths have equal temperatures. We begin by ex- 
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FIG. 3. Thermal radiation |£s_| 2 for the system described in Fig. [U normalized by the maximum of the 
internal and external bath temperatures T = max { T d . T c }, as a function of the dimensionless frequency 
du/'y where dcu = (to — coq), with fixed A = —2.5 and under various operating conditions. The radiation 
spectrum is shown (a) in the limit T d T e of negligible externally incident radiation, for different values 
of C and (b) for fixed £ = 4 but different T d . T e and linear decay rates 7d, 7e- 


ploiting a simple perturbation theory in which the thermal fluctuations of the cavity-field da and 
radiation Ss- = — £ e + ^2^5a are treated perturbatively, leading to analytical expressions for the 
corresponding thermal energy and radiation spectra. Assuming |s p | 2 7 A: B {X d . T e }, we obtain: 

/IX 1 2 \ & B /+(27 d T d + 27 e T e ) 

(H ) =-- 

<|fc_| 2 } = k B T c + i^[ 7d (/ + T d - f_T.) + 7c (/+ - f-)T e \ 

where / + , /_ and D are given by 

f± = (u + u 0 - 2o/p - 2ct|a 0 | 2 ) 2 + 7 2 ± a 2 |a 0 | 4 
D = [7 2 + (cu 0 - w p — 2a|a 0 | 2 ) 2 — (u> — a/ p ) 2 - a 2 |a 0 | 4 ] 2 
+ 4 7 2 (a :- u v )\ 

and where |a 0 | 2 denotes the steady-state cavity energy in the absence of fluctuations, the solution 
of ©. Here, for simplicity we have assumed the classical limit @(cuo, T) —> k\>T. In the absence 
of the external bath T e = 0, the above equations are similar to those obtained in the case of me¬ 
chanical oscillators.— The situation changes with the channel due to the unavoidable intereference 
and induced correlations of the emitted and reflected radiation, described in (7]) by the / terms. 


( 6 ) 

(7) 
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Figure [3] illustrates the radiation spectrum of the cavity under different operating conditions, 
showing excellent agreement between the numerically (noisy) and analytically (dashed lines) com¬ 
puted spectra. We note that all of the results shown in Fig. [3] correspond to cavities operating out¬ 
side of the bistable regime: although it is possible to obtain a complete and analytical description 
of the spectrum based on ([6]) and ([7]). such an analysis is difficult in the bistable regime due to 
hopping between states, requiring a complicated description of the transition rates and stationary 
distributions of the system.™ For instance, in the bistable regime, one observes the appearance 
of a temperature-dependent supemarrow spectral peak whose amplitude decreases with increas¬ 
ing £. In the absence of bistability, a similar but weaker amplification occurs as ( —>■ |A/2| or 
equivalently, as the cavity frequency becomes resonant with the pump. Regardless of regime, at 
suffiiently large ( (once the cavity resonance has crossed u; p ), modulation of the thermal noise by 
the pump causes the spectrum to transition from being singly to doubly resonant due to the emer¬ 
gence of an additional anti-Stokes peak.— In the limit as ( —> oo, both peaks move farther apart 
and their amplitudes assymptote to a system-dependent constant. 

Figure Ob) explores the dependence of the peak amplitudes on various cavity parameters, in¬ 
cluding T d = T e , T d T e , and T d -C T e , corresponding to a resonator that is either at thermal 
equilibrium, heated, or cooled with respect to its surroundings, respectively. When noise entering 
the system through the external bath is negligible T d T e , similar to the previously explored 
situation involving mechanical oscillators,— one finds that the Stokes peak is always much larger 
than the anti-Stokes peak (green line). Essentially, for a > 0 the cavity nonlinearity favors down- 
conversion, as captured by the asymmetric /+ terms above. The peak radiation associated with 
the Stokes peak can be readily obtained from © in the non-bistable regime A < — \/3, and is 
given by max |<5s_| 2 = ^^(1 + 2A 2 )k B T d , reaching 7k B T d precisely at the onset of bistability 
and when y e = 7d- At larger ( > | A/21, the amplitude of both peaks decreases with increasing ( 
where, as ( —> oo (not shown), the amplitude of the Stokes peak —)■ k B T d while the anti-Stokes 
peak dissapears. The situation changes dramatically when the noise entering the system through 
the external bath cannot be ignored, i.e. T e > T d . In particular, as observed from ©, although 
the cavity spectrum favors Stokes to anti-Stokes conversion regardless of the relative tempera¬ 
tures or decay rates, we find that the spectrum of outgoing radiation can be dramatically different 
depending on the regime of operation. When T e T d where noise is dominated by external ra¬ 
diation, we find that the anti-Stokes peak dominates (red line) except when y e v d , at which 
point the spectrum exhibits a symmetric lineshape (orange line). Such a reversal in relative am- 
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plitudes is captured by the /_ terms above, which include correlations and interference between 
the emitted and reflected radiation. The maximum radiation in the non-bistable regime in this 
case max |<5s_| 2 = k^T^l + Af-A 2 — reaching 25k B T e at the onset of bistability and when 
7e 3> 7d (rather than under rate matching). Interestingly, we find that when the two baths lie at 
the same temperature T e = T d , both peaks have equal amplitudes regardless of 7 e /7d, though the 
maximum amplitude in this regime also occurs in the limit y e This unexpected symmetriza- 

tion of the spectrum arising due to interference effects seems to be a unique property of thermal 
radiation in this system. Although previous work on nonlinear fluctuations in the quantum regime 
observed similar peaks in the spectrum, a symmetric spectrum was found to arise only at zero 
temperature (a singular point of the theory^) due to quantum tunneling.— 

Although a number of the abovementioned effects have been observed in mechanical oscil¬ 
lators, they remain unobserved in the context of thermal radiation where they could potentially 
be exploited in numerous applications: 2 - ' 25]26 As demonstrated above, the interplay between the 
internal and externally incident radiation and the coherent pump leads to new effects in thermal 
radiators, including dramatic changes in both the coherent and thermal output spectrum with tem¬ 
perature, along with temperature-tunable Stokes and anti-Stokes side peaks. Finally, we conclude 
by proposing a realistic, silicon ring resonator design, depicted schematically in Fig. Q] where 
one could potentially observe these effects near room temperature and with operating Q ~ 10 5 
and input power k P l 2 - lmW, leading to a\s p rQ 2 ~ |A|. For these parameters we find that 
IsTF Asm ^ ~ 0.04K -1 at T ^ 300K. Although this is almost two orders of magnitude smaller in 
comparison with thermo-optic effects in silicon, which lead to tunable powers ~ K” 1 for the same 
structure, at lower temperatures T < 100K where the thermo-optic coefficient is much smaller—, 
our fluctuation-induced effects offer significantly better temperature tunability. Other cavity de¬ 
signs such as the nanobeam cavity described in Ref. [35! yield much larger a and allow smaller Q to 
be employed, leading to even larger tunability compared to that obtained via thermo-optic effects. 

We are grateful to Mark Dykman for very helpful comments and suggestions. This work was 
supported in part by the National Science Foundation under Grant No. DMR-145483. 
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